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1 Introduction 

Advances in observational cosmology in the last decades have provided us with countless 
ways to test cosmological models and so far, the Standard Model of Cosmology has passed 
all of them. In the simplest version this standard model of cosmology comprises a period 
of accelerated expansion, inflation, at the very beginning and has the dynamics of the uni¬ 
verse dominated by dark matter, and later dark energy, after an initial epoch of radiation 
domination. 

In recent years the results from the WMAP [1] and the Planck [2] missions have pro¬ 
vided us with an unprecedented amount of information from the temperature anisotropies 
in the Cosmic Microwave Background (CMB), constraining the allowed parameter space. 
Also, recent B-mode polarisation measurements [3] have put additional tight bounds on the 
tensor-to-scalar ratio. In spite of all this progress, the correct model of the early LTniverse is 
still unknown, as there are many inflationary potentials that are consistent with the experi¬ 
mental results. Future surveys will provide an even further improvement in the experimental 
accuracy and should probe more scales than ever before. However, in order to make the most 
of these new observations, further efforts are required also on the theory side, for example to 
compute observables with the precision required by the experiments. 

According to inflationary theory, the seeds of structure were generated during inflation 
through the quantum fluctuations of one or more light scalar fields. To study the evolution 
of these initial perturbations we use Cosmological Perturbation Theory. This amounts to 
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adding inhomogeneous perturbations to the Friedmann-Lemaitre-Robertson-Walker (FLRW) 
solution and studying their evolution equations. While the linear theory seems to have been 
successful so far, higher order corrections seem necessary, in order to study new observables, 
such as uon-gaussianities, and to make more accurate predictions of the observables we 
already measure. 

The gauge invariant curvature perturbation on uniform density hypersurfaces, £, is a key 
variable when computing predictions from inflation [4-7]. Its correlation functions encode the 
relevant information about the micro-physics of inflation and are also related with observables 
we can probe, such as the temperature anisotropies in the CMB and the distribution of 
galaxies in the sky. It is, therefore, fundamental to have a precise knowledge about the 
evolution of this curvature perturbation. It is well known that, on super-horizon scales and 
for adiabatic perturbations, C, is conserved at first order (see e.g. Ref. [8]). This result has 
been extended to higher orders using non-perturbative methods such as the 5N formalism 
[9] and using other approaches (see e.g. Refs. [10, 11]). At second order in cosmological 
perturbation theory, the conservation of the curvature perturbation has been established for 
scalar perturbations in Ref. [12]. 

However, the evolution equation at all scales and at second order in the perturbative ex¬ 
pansion including scalar, vector and tensor perturbations has so far not been fully explored 1 . 
Furthermore, at second order, there exist different definitions of this gauge invariant per¬ 
turbation, and it is not clear that the conservation of one implies the conservation of the 
others. 

In this work, we aim to address both of these issues. We start by reviewing the different 
versions of the gauge invariant curvature perturbation on uniform density hypersurfaces and 
show how they are related. We then derive the evolution equation for each convention and 
compare the results. Besides the scalar contributions we also keep all vector and tensor con¬ 
tributions, as well as the anisotropic stress. Finally, we take the large scale limit and check 
for the conditions of existence of conserved quantities. 

The paper is structured as follows. In the next section, we present the different con¬ 
ventions for the metric perturbations and give the necessary gauge transformations. The 
different definitions of ( are given in Section 3, along with a number of auxiliary gauge in¬ 
variant quantities. A derivation of the evolution of (2 is presented in Section 4. We then 
present our conclusions in Section 5. 

2 Definitions and Conventions 

In this initial section, we introduce our conventions for the metric and stress-energy tensors. 
We simultaneously perform the 3+1 decomposition of space-time and and also the scalar- 
vector-tensor (SVT) decomposition, in which tensors in the spatial slicing are split into their 
scalar, divergence-free vector and transverse and traceless tensor parts according to their 
transformation behaviour on the 3-dimensional submanifolds. For the most part, we follow 
the notation of [17]. All quantities T are expanded as 

T = T 0 + ST, + Ut 2 + ..., (2.1) 

1 Note, however, the early works of Tomita in Refs. [13-15] regarding the interaction between scalar, vector 
and tensor modes in the Einstein-de Sitter universe, as well as more recent work in Ref. [16] and references 
therein. 
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in which the subscript denotes the order in perturbation theory. We will drop the subscript 
for the background quantities as the unperturbed ones only appear at their definition. 

Throughout the paper we assume a FLRW background spacetime with zero spatial 
curvature, and use conformal time. Greek indices, /jl, v, A, range, from 0,... 3, while lower 
case Latin indices, i,j, k, denote spatial indices ranging from 1,... 3. 

2.1 Metric Tensor 

The metric tensor can be split in several different ways. We will show four distinct ones, 
which vary in the way the spatial part of the metric is arranged. The version that we will 
use in most of the calculations below takes the following form, at all orders 

g 00 = ~ a * 2 (1 + 2(f)) , (2.2) 

gto =a 2 (B ti - Si) , (2.3) 

9ij =a ~ J [$ij + 2 Cij] , (2-4) 

in which cf> is the perturbation to the lapse, B and Si are, respectively, the scalar and vector 
parts of the shift and Cij is the perturbation to the spatial part of the metric. The first 
convention we will treat is defined by arranging Cij as 2 

Cij — T C.ij T T hij , (2.5) 

in which if} is the curvature perturbation in this metric convention (as used by Mukhanov, 
Feldman and Brandenberger in Ref. [18] and Malik and Wands in Ref. [17], for example), E 
and Fi are, respectively, a scalar and a vector part of the spatial metric and hij is the tensor 
potential, representing gravitational waves. The vector and tensor quantities above obey 

S\ t = 0, F\ i = 0, h 13 ■ = 0, h\ = 0. (2.6) 

This first convention for ip can be understood, at first order, as the perturbation to 
the intrinsic curvature, as explained in Appendix C. The other conventions do not have this 
property, but can be generally understood as perturbations to the scale factor aft). Appendix 
C also contains a definition of the scale factor from the extrinsic curvature, which is more 
easily relatable to the versions of vp given below. 

A variation from the form (2.5) consists of collecting the trace of Cij in a single variable, 
here denoted by ipT- This split was used, e.g., by Bardeen in Ref. [19] and also by Kodama 
and Sasaki in Ref. [20], where V’T was denoted by H^. The perturbation to the spatial part 
of the metric becomes 

Cij = —ipT&ij + Bjj — —dij'V^E + F)jj) + hij , (2-7) 

which, upon comparison with the previous convention, (2.5), shows that the new curvature 
perturbation is given by, at all orders 

ipT = - ^V 2 F. ( 2 - 8 ) 

2 Note that the definition of the tensor perturbation is slightly different from the one in Ref. [17], i.e. hC lw ' > = 

2 hij . 
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The third kind of decomposition of gtj we will treat is similar to the second one, Eq. 
(2.7), but factors out the determinant of the spatial part of the metric, instead of the trace. 
This is the decomposition used by Salopek and Bond in Ref. [5] and also by Maldacena in 
Ref. [21]. It can be written as 


g lj =a 2 e^ D le“] i j 


(2.9) 


in which w is a traceless tensor and is the curvature perturbation of interest in this 
convention, defined by e 6 ^ D = det (gij/a 2 ). This quantity is usually interpreted as being 
a perturbation to the number of e-folds [5], N, given by IV = In a — V’zn A related 

interpretation would be to think of it as a perturbation to the volume of spatial hyper surfaces, 
as it is proportional to the determinant of the spatial metric. It can be shown [17], that the 
first and second order parts of iJjjj are related by the following expressions, in the notation 
of this work 

VtDi = - Vti = -Vd + ^v 2 Ei , (2.10) 

i’D2 = — i’T2 — T)p\ijC l \ = (2.11) 

= - v>2 + \v 2 E 2 - 2ipl - | h U jti{ + ^iV 2 Ei 

- f ft? + Fuj) - ( 2 f; j + Ef) . (2.12) 

The fourth convention is not a variation of gij per se, but only a different way of 

defining the curvature perturbation. As with the third definition, Eq. (2.9), we factor out 
the determinant of the spatial part of the metric, but in this case, we use the inverse metric 
to do so. Therefore, it is now defined as 


g ij = a- 2 e- 2,/ "[e^P, 


(2.13) 


in which, again, up is a traceless tensor and i>j is the new version of the curvature per¬ 
turbation, determined by e -6 ^’ 7 = det (g^a 2 ). To our knowledge, this is the first time this 
definition has been used in the literature. Concerning its interpretation, it can still be seen 
as a perturbation to the scale factor and we find it to be equal to the integrated expansion, 
when the latter is evaluated in a comoving threading (see Appendix C for more details). 
Comparing this new version of i> to the original one, we find the following relations 

i’ll =~ipi + ^V 2 £q, ( 2 - 14 ) 

'ii2 = ~ ^2 + ^V 2 £ 2 - 2 ii 2 - h Uj K j + \i>i V 2 Ed + ^ (B lti - S u ) (B( - S*) 

- (Ei,H + Fuj) - §*1 - \Ei tij (2 Ft* + Ef) . ( 2 . 15 ) 

We will use these four conventions to define different versions of the gauge invariant 
curvature perturbation in the next section. 


2.2 Stress-energy Tensor 

As for the stress-energy tensor, it is defined in the so-called energy frame [22, 23] 3 , including 
anisotropic stress: 

T/j,,, = (P + p)u tl u u + Pg^ + 7 ^ , (2.16) 

3 This frame is defined by the condition that the 4-velocity, i/ 1 , is an eigenvector of Twith eigenvalue p. 

This is equivalent to setting the energy flux, to zero. 
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with 


7r /J M = 0 , and tt = 0 . (2-17) 

The perturbative expansion is the standard one, as given in [17]. We merely reproduce the 
SVT splitting of n 2 ^ as it is non-trivial at this order: 


71200 = 0, 7T2iO = ~^1 ij (vyi + V ’l) > 


^"2 ij — ^ 


II2 ij + n 2(ij) + Il^jj 


- -<%v 2 n 2 


+ — 7T ifc/ C*! , 


(2.18) 


where the quantities II, 11* and II l? are, respectivelly, the scalar, vector and tensor parts of the 
linear piece of the anisotropic stress tensor, v and Vy are the scalar and vector components 
of the velocity perturbation and C\ij is the spatial metric perturbation, as given in equa¬ 
tion (2.5). The terms that arise in 7T2 *o and 7 T^ij which depend on the first order anisotropic 
stress, 7Tiij, are essential in order to not violate the constraints that define the anisotropic 
stress, Eqs. (2.17). Should these terms not be included, second order contributions of the 
anisotropic stress to the evolution equations are not well taken into account. That seems to 
be the case for the energy conservation equation given in [17], which differs from our result 
presented below in equation (4.4). 


2.3 The gauge problem 

Most of the quantities defined above will depend on a choice of coordinates on the perturbed 
manifold, the so-called gauge choice [17-20, 24]. A gauge generator vector field £ is defined 
to parametrise the change of coordinate system. The gauge transformation it induces in a 
tensor T is given by 

T = e £ «T, (2.19) 


with £^ being the Lie derivative in the direction of £. The transformation rules for the 
perturbations are obtained by expanding the relation above order by order. Up to second 


order, the gauge transformations are given by 

<m=5Ti + %T 0 , (2.20) 

6 T 2 = ST -2 + %T 0 + £\T {) + 2%<m , (2.21) 

where the gauge generator was also expanded order by order as + ■ ■ • • It can 

be decomposed further into scalar and vector parts as 

(n = M’*+y)- ( 2 . 22 ) 

Applying these transformations to the metric tensor, one finds the first order metric potentials 
of the original convention, Eq. (2.5), to change to 

4>i = (j>i + £Lai + a[ , V’l = V’l _ % Q i, (2.23) 

E\ = E\ + Pi, B\ = B\ — a\ + (3[ , (2.24) 

Fi = Fi + 7 *, 5{ = 5i- 7 f, (2.25) 

h\ j = ]£( , (2.26) 
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and the first order fluid quantities to transform to 


(2.27) 

(2.28) 


dpi = dpi + a\p , 8 Pi = 8Pi + a.iP ', 

vi = vi - /3[ , v l vl = v l vl - a. 

The gauge transformations at second order are given in Appendix A. 

The ambiguities created by this gauge issue are solved by “fixing” or choosing a gauge 
and hence constructing gauge invariant quantities [12, 17-20]. In the following section we 
derive gauge invariant quantities, specifying a particular gauge by choosing convenient tem¬ 
poral and spatial hypersurfaces, or specifying the slicing and threading, respectively (for a 
non-technical introduction to the gauge issue in cosmology see Ref. [25]). The gauge invari¬ 
ant quantities are uniquely defined by specifying the hypersurface (at all orders), and can 
therefore also be evaluated in other gauges, if necessary. 


3 Gauge invariant quantities 


The method we use to generate gauge invariant variables starts with performing a gauge 
transformation on a variable of interest, e.g. tpi. One then substitutes the gauge generator 
components Cf with those obtained by solving a gauge fixing constraint, e.g. 8 p 1 = 0. The 
end result is a gauge invariant quantity, e.g. the curvature perturbation in uniform density 
hyper surfaces, Cl- We apply this method for the quantities of interest in the subsections 
below (for details see e.g. Refs. [17, 24]). 


3.1 Curvature perturbation on uniform density hypersurfaces 

The focus of this work is the curvature perturbation on uniform density hypersurfaces (. As 
was already sketched above, it is defined to be equal to —ip in the gauge in which the density 
field is uniform ( 8 p = 0). Starting with our first convention for the metric, Eq. (2.5), this 
condition is sufficient to fully construct ( at first order as ([4, 8]) 

Cl = (3.1) 

P 

However, at second order, one is also forced to specify the first order gauge to define this 
curvature perturbation unambiguously. For this convention of the metric tensor, Eq. (2.5), 
we will use the following gauge conditions to define C 2 ([12]) 

C 2 = -^2 , if 8 p 2 = 8 pi = Ei = 0 , F[ = 0 . (3.2) 


These add a flat threading to the uniform density gauge. The general formal expression for 
C 2 is given in [17]. In full detail, the formula is rather complicated and we write it here with 
the r.h.s. evaluated in flat gauge 4 , 


H 


P 


2 U 


,k 


C 2 — j 8 p 2 4 -pi ( 2 R + T-L' — H— ) 8Pi 4 j^^Pi^Pi — iry^^Pi^Pi 7 (-®i ,k ~ Sifc) 8 p\ 


P- 


+ V 


-2 


,Rch) 


^2 8 p'l 8 pi + (Spi'Bf - 8 pfS[ 


2 p'- 
1 


P 


4 7 [h-i + 2 Phi J 8p X i 

8j P 


4 Flat gauge is defined by the conditions tp = E = F l = 0. 
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We can see that, in contrast with the first order result, the second order ( is much harder to 
relate to density perturbations in flat gauge, given the presence of vectors and tensors. In 
spite of this, this expression is still useful in writing the gauge invariant curvature perturba¬ 
tion in terms of multiple scalar fields, as is done in [26, 27]. 


Let us now move to the second convention of the metric, Eq. (2.7). In this case, 5p = 0 
is no longer a sufficient gauge condition to define an invariant, even at first order; one must 
also specify the scalar part of the threading, due to the inclusion of E in the definition of i/jt 
(see Eq. (2.8)). The extra condition that is most often chosen is v\ = 0, which results in the 
following expression 0 

Cti = -V’Ti - + ^V 2 [ vidr , (3.4) 

P 3 J 

in which the integral in conformal time is indefinite. The introduction of these integrals is 
the disadvantage of using the gauge condition, v\ = 0. This might be problematic, as this 
condition only sets the gauge up to an arbitrary function of the spatial coordinates, which, 
in turn, might spoil the gauge invariance of the new variable. In spite of this, it is possible 
to construct a gauge invariant quantity, by defining it to be 


Cti = Ci + ^V 2 / Jidr . 


7 ' 


with J\ being the gauge invariant velocity on flat hyper surfaces, defined by 

J\ = Ei + v±. 


(3.5) 


(3.6) 


While the integral in Eq. (3.5) is still indefinite, the integrand is gauge invariant and, there¬ 
fore, this is the definition we use. 

At second order, one sets the second order gauge in the same way, i.e. bp 2 = v 2 = 0 
and, to avoid additional issues with indefinite integrals, one can choose 5pi = E\ = Ff = 0 
for the first order gauge fixing. With this choice, we find 

Ct 2 = C 2 + V J J 2 dr , (3.7) 

in which J 2 is the second order equivalent of J\ in this gauge, i.e. it equals E 2 + v 2 in the 

gauge obeying 5pi = E\ = F[ = 0. As is visible in the expression above, Eq. (3.5), the 
only variable of interest is V 2 J 2 and hence, for shortness of presentation, that is all we show 
below, with the r.h.s. evaluated in flat gauge 


2 r 


Li 


V J 2 —V V 2 + — 5Pi \ T~L{yyi + v'l ) — Vyi — v{ 


+ 


2 P- 


/2 


+ 


&PlM , 


(3.8) 


+v 


-2 


■^ 2 6 Pi S Pi + ^ ($P'i B 1 ] - tpfsf) - ^ (hf + 2Uh l {) Sp hij 


- 1 


As we will see in Section 4, this quantity is relevant regardless of the choice of convention for 
the metric, as it will appear in the evolution equation for the curvature perturbation. 


5 An alternative ch oice would be E\ = 0, but that would simply result in the expression for the original 

metric convention, as i/>ti = V ’ l - 
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Let us now turn to the third convention of the metric, Eq. (2.9). For this case, (r> will 
be defined as being equal to V’D instead of — i/jp, in order to keep the same sign as £. Starting 
at first order, we see that we get either Cdi = Cl or Cdi = Cti, depending on whether we 
choose E\ = 0 or v\ = 0, respectively, for fixing the threading. The second order result 
is more interesting, as there is no gauge fixing for which it is equal to either of the other 
definitions above. In the most conservative case, the choice of gauge fixing is bp 2 = E 2 = 0 
at second order and bp 1 = E\ = F\ = 0 at first order. This results in 1 ' 

Cd2 = C 2 - \h U] h\ 3 - 2(1 . (3.9) 

A different gauge fixing is 5p 2 = = 0 and bp l = v\ = vyi = 0, for which the result is 

Cm =C 2 + \J V 2 J 2 dr - h Hj h\ 3 - 2(f (3.10) 

+ 2 Ci,i J ( J( + VI) dr + | J | [( J? + Vl) Ti] . + V 2 J ( J’i + Vi) dr' j dr , 

in which V[ is the gauge invariant velocity vector perturbation in flat hypersurfaces and T is 
the gauge invariant lapse perturbation in uniform density hypersurfaces. In a general gauge, 
these quantities are given by 

Vi = V vl + Ff, (3.11) 

Tl (3.12) 

For the forth version of the curvature perturbation, Eq. (2.13), the procedure is very 
similar to the one for the third convention. As in the previous case, the first order quantities 
obey (i i = Ci or ( n = Crii depending on whether E\ = 0 or v\ = 0 is chosen for setting the 
threading. At second order, the results are 

C /2 = C 2 - ~ 2Ci + l (W u - Fh + A hi - J M ) (W{ - Vi + A( - jfj , (3.13) 

if the gauge is fixed with 5p 2 = E 2 = 0 and 5p x = E\ = F\ = 0, and 

=c 2 + U v 2 J 2 dr - hi H] h\ 3 - 2 C 2 + ^ (Wu + (wi + Afj (3.14) 

+ 2 Ci,i J {J( + Vi) dr + ^j | [( Ji + Vi) Ti] . + V 2 J u J {J( + Vf) dr'j dr , 

when the gauge choice is bp 2 = v 2 = 0 and 5pi = v\ = vyi = 0. The new first order 
gauge invariant quantities that appear are the vector velocity in zero shift gauge, W{, and 
the momentum perturbation in uniform density gauge, A\. They are given by 

W[ = v i y 1 -S[, (3.15) 

Ai=ui + Bi + %. (3.16) 

P 

6 This result is well known in the case without tensors. See, for example, Refs. [7, 11, 27]. 
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3.2 Non-adiabatic pressure 

One of the quantities determining the evolution of the curvature perturbation is the non- 
adiabatic pressure [8, 19, 20, 28, 29]. It is defined as the deviation from the adiabatic relation 
as 


5P = c 2 s 5p + 5P nad , (3.17) 

with Cg the adiabatic sound speed defined as c 2 = P' / p '. At first order, this definition 
automatically generates a gauge invariant quantity, but, at second order, this is not sufficient 
and many choices can be made. Our first choice is to define it as the gauge invariant quantity 
that reduces to Eq. (3.17) in the gauge in which 5p 1 = E\ = F[ = 0. In a general gauge, this 
quantity is given by 

2 f p" p' 0 "\ 

8Pna d2 = SP r 2 ~ Cg^P2 — ~i8p\8P\ + y~y2 - 

+ ^-SpM - 2 (Fi + E’f) AP nad i,i. (3.18) 

With the different choice of threading, v\ = vy{‘ = 0, one finds instead 

SP^ 2 = SP nad 2 + 2 5P nad i,i J (yi + jf) dr . (3.19) 

For a barotropic fluid, with P = P(p), both expressions vanish, as can be easily checked by 
evaluating them in their defining gauge, i.e. with 6 p 1 = E\ = F{ = 0. 

The quantities presented so far include the full set of gauge invariant quantities required 
for the full derivation of the evolution equations below. 

4 Evolution equations 

In this section, we present the derivation of the evolution equations for all versions of £. Our 
strategy consists of calculating the derivative of expression (3.3) and using only the perturbed 
energy-momentum conservation equations up to second order to simplify the result. Lastly, 
we substitute the gauge dependent variables for gauge invariant ones, using the expressions 
found in the previous Section, to arrive at our final result. Having found the result for (2 
in the original convention of the metric, Eq. (2.5), we then rewrite the evolution equation in 
terms of the different definitions of (. 

4.1 Fluid equations 

Energy-momentum conservation, T 1 '^ = 0, governs the evolution of the fluid density and 
velocity. For simplicity, we present these evolution equations in flat gauge. The first order 
energy conservation equation is given by 

8p[ + m (8 Pl + <ypi) + (p + p) vV = 0 , ( 4 . 1 ) 

while momentum conservation is 

8 Pi,k + (p + P) [Z[j t + 01 ,k + (l — 3c^) PZik] + —V 2 ni ) fe + — V 2 nu, = 0 , (4.2) 
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where the momentum perturbation Z is given by 

Z k = v v f - S k + B* + v{ k . 

At second order, we only require the energy conservation equation, which is 

Sp' 2 = -m ( 5p 2 + 5P 2 ) -(p + P) V 2 u 2 - 2 (<5Pi + 5/q) vV - 25p ltk (v^ + v\ k 

- 26P lik Z 1 k -(p + P) [±Z' lk Z k + 2(1- 3c 2 ) UZ lk Z k + 2 ^ k Z k 

+20i , k (wf + + ,A: ) + 20 1 v 2 u 1 - 4 /i' Uj v] - Zj* Qv 2 n ljfc + v 2 n lfc 

- 2 (ny + + u/ j - i^v 2 nx 


(4.3) 


(4.4) 


The above equations are sufficient to derive evolution equations for the curvature perturbation 
at first and at second order [ 8 ]. 


4.2 Evolution of the curvature perturbation 

We can now derive the evolution equation for the curvature perturbation on uniform density 
hypersurfaces. We follow the strategy stated at the beginning of this section. At first order, 
the result is well known to be 

Ci = 4 v +-*yfff > ( 4 - 5 > 

where only the first order energy conservation equation was used. On large scales (“V —» 0”) 
and in the absence of non-adiabatic pressure, one finds the familiar conservation equation 

c; = o. 

Ay) Ay) 

For the other conventions for the curvature perturbation, + D [ and Qj , the evolution 
equation at first order is the same and is given by 


Cti 


-H 


&Pnad 1 
P + P 


(4.6) 


which shows these versions of (4 are conserved at all scales, when non-adiabatic pressure is 
negligible [9, 30]. 

At second order, the complexity increases. The procedure to obtain the final result is as 
follows: use the energy conservation equation at first (Eq. (4.1)) and second order (Eq. (4.4)) 
to substitute for Sp\ and bp' 2 and substitute (IV 2 !!^ + V 2 ]!^) with the momentum con¬ 
servation equation, Eq. (4.2). The last step is to use the defining expressions of the gauge 
invariants to eliminate all gauge dependant variables. The final result is given by' 


— C2 + 2£i — — iWu + +!,«) ( W\ + A) J + ] — 


(4.7) 


1 . 


-y z j 2 + u 


2 n 


^-P na d 2 

p + p 


-2 n 


<TLu . ^ 

(,P + P ) 2 3 

1 


T ( VI + ) j . + 2Ci,i ( VI + 


- — (n Uj + n lM + u Uj - -b ij v z n,) [h? + v ^ + j? 


"Note the absence of inverse Laplacians. That is explained by an exact cancellation between the terms in 
(2 and those in V 2 J 2 , as can be shown by comparing equations (3.8) and (3.3). 
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We are now able to identify the different terms that source the evolution of ( 2 - We note, in 
particular, the appearance of vector and tensor source terms as well as the anisotropic stress 
which did not appear at first order in this equation 8 . 

We are now in the position to substitute for the other versions of C and find their 
evolution equations. For ( 2 T, we find 


(—C 2 T + 2Cf — - [W\i + A\^ (w[ + + -hu-jh'\ \ — (4.9) 


n ; ^-Pnad 2 _ ^ . ^nad 1 , ^ 

P + P H (p + P) 2 + 3 


r^ + Jf)] + 2 Ci,i (vi + ,/f) 


y (niy + n 1(iii) + n Mi - ^v 2 n x ) (hf + vf* + j ( j ), 


while Cd 2 evolves as 


Cd 2 ~ l (Wu + Ai,i) (W[ + A{ 


-v 2 j 2 + n 6Pnacl2 _ 2 n 6P ™ A1 


p + p 


{p + P) 2 3 L 

i 


- y (nib + n Kid) + n bb - ^bV 2 n^ (h\ jl + vi’ j + J( j ) , 


and the result for Q D2 is 


-+ - \ (u'l. + d,,) (w';-.4/)) 




M §P 2 


nad 2 _ 2T~L 


nad 1 


Art ( 1 

Y f nr 1 iy + + ni p — -5 tJ v rii 


p + p ““'(p + P) 2 

i +vi' j +j( j ). 


(4.10) 


+ — 4" (Vi + J\ ^ + 2 Cdi,i (V\ + J[ ^ 


(4.11) 


(i;) . 

The simplest evolutions equations are found for the C /2 and Q 2 versions of the gauge invariant 
curvature perturbation. They are given by 


-C /2 + 77 vi i + J\,)(v; + j?- 2 w{- 2 a; 


.* - 


(4.12) 


3 p + P 


- 2 n 


^nadi 2 
{p + P) 2 3 L 


r (V{ + j( ) +2C nAvi + Ji 


C 1T~L ( 1 

y ^niij + + iii ,ij — rii 


+ vl' 3 + J 


,v 


8 Note however, that the scalar part of the anisotropic stress tensor would source the evolution of C at first 
order by acting on the evolution of V 2 J. This can be seen more clearly by deriving Eq. (4.5) and using the 
momentum conservation equation, Eq. (4.2), to substitute for V 2 J: 


Ci+nCi- 


p’ 


3 (p + p) 


V 2 Ti-iv 2 4>i + 


H 


4 T|i;i/I 1 
p + p 


+n 


; 4P na dl V 4P na dl 2 
p+P 3 (p + P) 9 (p + P) 


V 2 V 2 ni =0, (4.8) 


in which <f>i is one of the Bardeen potentials, given in terms of our variables as lii = Ti+'H(Ti —Ji) + (2li —Ji)'. 
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and 


A v )r 

V 2 

zu 

p' 


xp(' u ) Xp2 

_ •p__nad_2 ^Wiadl 


P + P 


(P + P) 2 

+ III ,ij ~ V 2 IIi ^ 



(4.13) 


This final expression, like its first order version, Eq. (4.6), shows that, in the absence of 
non-adiabatic pressure and anisotropic stress, this version of the curvature perturbation is 
conserved on all scales. While this is interesting, in order for this result to be useful, one 
would likely be forced to estimate the integrals in the defining expression for (j 2 , Eq. (3.14). 
This is not likely to be straightforward, given the indeterminate nature of the integrals. This 
evolution equation matches the results of Ref. [31] for the integrated expansion in the absence 
of anisotropic stress, obtained in the covariant approach. 


4.3 Large scale approximation 


Here we perform the large scale approximation, by neglecting all terms with spatial derivatives 
in the equations above 9 . This limit simplifies the evolution equations to 


1 


-C 2 + 2Ci - ^w u wi + -h Uj h l { = n 


,$Pi iad_2 

P+ P 


SP'l 


nad 1 


-CD2-\w u wi] =n 5 -^^--2 n 

3 ) P + P 

(v) 

for the evolution of both C,D 2 and and 

-Ci 2 + \v u (Vl-2Wl)) =n d -^^-2n 

3 J p + P 


(v)/ _ ^w ( l-Pnad2 


-C " = n 


P + P 


-2 n 


~"'(P+P) 2 

and 


^nadl 

m T 

(P + P) 2 

~ P' 

, <„di 

2H 

'( P+P) 2 

pf 

r ^nadi 

2 n. 


(,p + p) 2 p’ 


(4.14) 

fif, (4.15) 

i(4-16) 
,hf- (4-17) 


Note that, in all cases above, the pairs are equal in the large scale approximation, except 
for (j2 and which have a different contribution from vector perturbations. We give the 
expressions for each version in this approximation in Appendix B. From this result, one can 
see that, even in the absence of the scalar non-adiabatic pressure, AP na d, neither curvature 
perturbation is conserved, 

(-C'2 + \hujtil - IwlWu) = (-CD2 - -W\W X ^ = -eg' = -yn Uj hf . (4.18) 

(u) 

However, if the traceless, transverse part of the anisotropic stress, Tluj, is negligible, Q 2 is 
in fact conserved 

eg' = (eD2 + \wiw u ) = (c 2 + \wiWn - = 0 • (4-19) 

9 This is generally well motivated in the case of some metric potentials, as one expects the perturbed metric 
to approach the background metric on large scales [9], and we will assume the same is true for the matter 
variables, including the anisotropic stress. Should this assumption not hold for the particular model under 
study, then the results in this section are not valid and one should use the full results from section 4.2. 
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Although is exactly conserved, the difference between C^' and C'd 2 on ly depends on 
vector perturbations, which are usually negligible. Moreover, using the vector part of the 
momentum conservation equation, Eq. (4.2), in the absence of anisotropic stress, we find the 
evolution of W\ is given by 

W' u + 1 - 3c 2 )W u = 0. (4.20) 

Thus, this vector perturbation is conserved during radiation domination (c 2 = 1/3) and, as 
a consequence, Cd 2 is exactly conserved during that epoch. In the general case, we may 
therefore write the evolution of Cd 2 on large scales as 

C,' m = -\u(l-3c 2 s )WiW u , (4.21) 

showing again that it may only have an appreciable evolution if the vector modes are large. 

The evolution equations simplify further in Einstein gravity, as, in the absence of 
anisotropic stress, tensor modes stop evolving and hence this new conservation law con¬ 
verges fairly quickly to the conservation of £2 itself. Hence, for Einstein gravity, all versions 
of the curvature perturbation are conserved up to second order on large scales, if both the 
non-adiabatic pressure and the anisotropic stress are negligible. However, should the evolu- 

(v) 

tion of vectors and tensors be appreciable, the version of £ which is conserved is Q , i.e., 
the version defined by the determinant of and by using a comoving threading to fix the 
gauge. 

5 Conclusion 

We obtained the evolution equation for the curvature perturbation at second order in cos¬ 
mological perturbation theory, valid on all scales. With the inclusion of vectors, tensors and 
anisotropic stress, this result allows for high precision calculations of correlation functions 
on all scales. We derive this for six different definitions of (, based on several different splits 
of the spatial metric and on various choices of the defining gauge. The results for the evolu¬ 
tion equations show a substantial difference in apparent complexity, being simpler when the 
threading defining ( was chosen to be the comoving one, i.e. v l = 0. Eq. (4.13) for the evolu- 
tion of Q 2 is particularly short, but its usefulness is unclear due to the existence of indefinite 
time integrals in the definitions of Q 2 and oP^ d2 - On the other hand, for the versions of 
( for which the threading was chosen with E = 0, or the original £ 2 , the definitions include 
inverse Laplacians (see Eq. (3.3)). In both cases, non-locality is present in some form, either 
in time or in space, and there is no version of the curvature perturbation which evades both of 
these issues. However, in both cases, the difficulties of the calculation are resolved by solving 
additional differential equations, both of which require boundary conditions. In the case of 
the inverse Laplacian, the equation to solve is a Poisson equation, which only depends on first 
order quantities at a single time, while for the case of the integrals in time, knowledge of the 
full time evolution of second order quantities is required (V 2 J 2 in Eq. (3.14), for example). 
This seems to render the quantities without integrals in time more amenable for situations 
that require the calculation of C, from its definition, such as when its value is evaluated from 
the value of scalar field or density perturbations. In any case, all these issues disappear in 
the large scale approximation, for which the inverse Laplacian term in question has a well 
defined limit and the integrals vanish. 

Moreover, we found that, on large scales, the evolution of £ is sourced by the trans¬ 
verse traceless part of the anisotropic stress tensor, as well as non-adiabatic pressure. Both 
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quantities must therefore be negligible for any version of ( to be conserved. Furthermore, 
the version of the curvature perturbation which is exactly conserved is the one based on the 
determinant of g d and comoving threading, c!f'\ Eq. (3.14). Other definitions may evolve 
with the evolution of tensor and vector modes, should such an evolution be allowed by the 
theory of gravitation under study. For General Relativity, however, vector perturbations are 
usually very small and the evolution of tensor modes is negligible in the absence of anisotropic 
stress; therefore all versions of the curvature perturbation are approximately conserved on 
large scales. 

The results presented here are valid as long as the energy and momentum conservation 
equations, Eqs. (4.1), (4.2) and (4.4), are satisfied. This will be true if the stress-energy tensor 
is covariantly conserved, i.e. = 0, and the connection is the Levi-Civita connection 

(i.e. no torsion is present). This is the case in GR, but also in other theories, such as Massive 
Gravity and Bigravity [32, 33]. The latter theories are interesting in this context, as the 
tensor modes evolve differently due to the non-zero mass of the graviton [34] and therefore, 
Cd and Ci would be the only versions of the curvature perturbation that are conserved. 

Furthermore, the usefulness of these results may be extended to theories of gravity for 
which 0. This is possible if one can perform a conformal transformation to the 

Einstein frame and apply the same ideas to the effective stress-energy tensor that arises as 
the r.h.s. of the new field equations. The difference between our standard scenario and a 
modified one is that the effective matter quantities thus defined, would not have the same 
physical significance as the ones we use in this work. Therefore, in those modified situations 
it may be less trivial to clearly say when the curvature perturbation is conserved, as, e.g. the 
effective 6P n a d may not be negligible when the true matter perturbations are adiabatic. The 
same could apply to the anisotropic stress. 

Previous results on the subject of conserved quantities have not included anisotropic 
stress [31] and have either done the calculations fully in the large scale approximation [9] or 
used a different quantity [11, 35, 36]. 
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A Second order gauge transformations 

Here we present the second order gauge transformations for the quantities used in this paper, 
given by Eq. (2.21) when applied to the metric and the stress energy tensor. We show only 
the transformations for the variables defined in the original convention of the metric, Eq. 
(2.5), as the others can be easily obtained by combining the rules given here with the first 
order ones in the main text. 
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(A.l) 

(A.2) 


The transformations for the metric quantities are given by 

02 = V >2 - Ua 2 - \X\ + |V 2 X ij dj , 

E 2 = E 2 + p 2 + 4V --’V 2 X'-' ; j - \\7- 2 X\ , 


with Xij given by (see Ref. [17] for more details) 


Xij = 2 


77^ H-) o 2 + 77 (o ioq + 


+ 4 


a i ipiij X 27-LCuj) + Cuj,kt,i + Cukii ,j + X'lkjCi' ,i + 2 (Buaij + B\jOi^) 

+ 477aq (£hj + £ij,i) — 2oi,iOi,j + 2£ifc ) j£i fc j + a\ (Cuj + + (£i i,jk + £i j,ik) £i k 

+ £li,klil k J + Clj.fcC/ ,i + Cli a l,j + £lj a l,i ■ (A.3) 


The transformations for the fluid quantities are 

<5/02 = <5/92 + u 2 p + ai (p"ai + po\ + 26p) + (25p + p'ai) k (j5’ k + 7 i) , (A.4) 

5P 2 = 5P 2 + a 2 P' + ai (P"a 1 + P’af + 2 5P') + (2 5P + P'ai) k (p* + 7l fc ) , (A.5) 

E 2 = v 2 -P' 2 + V~ 2 X v k k , (A.6) 


with 


= i'u (20i + of + 277ai) — aqCij 

— CiCn,fc + Ci'Cii.fc ~ 2ai + 77ui,:) + 2rJij^^ — 2v\^u^ ■ (A.7) 


B Large scale limit of the gauge invariant quantities 

In this appendix, we supply the expressions for the different versions of the curvature per¬ 
turbation in the large scale approximation, i.e. when all spatial gradients are taken to be 
negligible. 

The large scale limit of C 2 and Ct 2 is given by 

77 1 / o"\ . 277 

C 2 = Ct 2 = 7<5p 2 H p ( 277 2 + 77 ; — 77— ) 5p\ -\ -p^Spidp'i , (B.l) 

P P \ P J P 

(v) 

while that for C,D 2 and Cp D2 is 

Cd2 = Cr> 2 = --^bp 2 + ^2 ~ ^ ^7*1 + -^y<5/9] 5p\ — -h\tjh\ J , (B.2) 

and the limit of (72 and (^ 2 is 

C/2 = — —fip 2 + ^2 ' ~~ ^~p^) ~ -huj/t/ + -.S'ijS] , (B.3) 

cl? = ~^P 2 + y 2 (n' - 77<5 P 2 + ^<5piVi - . (B.4) 

These expressions agree with the <51V formalism, where comparison is possible (see 
Ref. [9]). 
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C On intrinsic and extrinsic curvature 


In this appendix, we aim to clarify the relation between the different definitions of ip and the 
perturbation to both the intrinsic and extrinsic curvature of hypersurfaces of constant time. 
We begin by looking at the intrinsic curvature scalar. It is given by 

( 3 ) i? = R + - K 2 + K> UJ K I1V , (C.l) 


in which R^ u and R are the 4D Ricci tensor and scalar, respectively, n^ is the unit normal to 
the hypersurface, K^ u is the extrinsic curvature and K is its trace. The latter are given by 

K ^ = ~\\ £ nlU , K = - , (C.2) 


with 7 the induced metric, given by 7 ^ = g^ u + n^riy. The normal, n M , is perpendicular to 
all vectors in the tangent space of the hypersurface. It is therefore often convenient to choose 
coordinates such that n* = 0. However, this specific coordinate choice means that the usual 
gauge transformation rules are not obeyed, and for this reason, we shall also compute these 
curvature scalars using the 4-velocity, u^, to define the spatial hypersurface. In the latter 
situation, we will denote quantities with the superscript (u). 

We now present the calculations of these quantities up to second order in cosmological 
perturbation theory. The intrinsic curvature scalar is found to be 


S {3) Rl = —V 2 -01 , 


(C.3) 


5 {3) R 2 = 


AXj'Z-L q /~i rris~ikn . K, rimn r> s-iti s~rm r 

4 V ^2 — 0U1 n -h DUi mn, u l ,k ~ ZLy l m, 


krrmn 


1 Qsvmn ( ri k I r<k ri k k 

' ^^1 ^lmn, k ' ^1 k,mn *-'1 ™k> n ^ 


1 mk,n 
k rik 


y lkn,m 


1 a I rik riJ n 1 ri 3 rim k rik rimn 

+ ^ l °1 1 ,n °1 jk,^l m, — °1 n,m° 1 ,k 


(C.4) 


where Cij is the perturbation to the spatial part of the metric. The relation between ip and 
the intrinsic curvature is clear at first order, as they are related linearly. This is the reason 
why the perturbation ip is called the curvature perturbation. However, this is only true for 
the original version of ip, as given by the definition (2.5), since all other definitions include a 
contribution from the metric potential E. In any case, at second order, this simple connection 
between the intrinsic curvature and ip is lost, as there is no simple relation between any of 
our definitions of the curvature perturbation and S^R 2 - 

Performing the same calculation using the 4-velocity to define the spatial hypersurface, 
one finds instead a connection to the curvature perturbation on comoving gauge, 1Z, since 
the first order result for ( 3 )_r( m ) is 10 

< 5 ( 3 )4 u) = Av 2 [ip! - H(v 1 + Hi)] = 4 v 2 ^i . (C.5) 

a z a z 

At second order, however, the result is no longer related to the second order comoving 
curvature perturbation 1Z 2 in a simple way, i.e. 5^R^ ^ £ j\7 2 1Z 2 . This can be seen by 

10 Note that this quantity has the expected gauge transformation properties, since, being 0 at the background 
level (because of the assumption of flatness), the Stewart-Walker lemma [40, 41] dictates it to be gauge invariant 
at first order. Notice that this does not happen in the calculation with n. 
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evaluating in comoving gauge (v = B = v l v = 0) and comparing it with Jj-V 2 ?/^- In 

this gauge, the intrinsic curvature is given by 

§ (3) r (u) = § (3) R2 + ^ w . + s i,j w .. 5 ( C . 6 ) 


where Wi and Wij are linear functions of the metric potentials. It is clear that this is not equal 
to 4V 2 ^2, as there are no further cancellations that would recover that result. Therefore, 
one must conclude that none of our definitions of has a straightforward interpretation as 
the perturbation to the intrinsic curvature at an order higher than first. 

Moving now to the scalar extrinsic curvature, we start by noting that it is proportional 
to the local expansion (or when choosing the velocity 4-vector to define the 

spatial hypersurface). It is well known that the integral of the expansion along world lines, 
with respect to proper time s, can be used to define a local scale factor [9, 11]. This integral 
is defined as 

a = U ^X ds = ~\j Kds ' (c.7) 

and the local scale factor is given by e a . This interpretation is further supported by the fact 
that, at the background level, one has a 1 = 71. At first order, one finds 

M = -E[). (C.8) 

This variable has some similarity with our definition of i/jt , but still has a contribution from 
B, which is not present in any of our versions of the curvature perturbation at first order. 
Turning now to the situation with u as the normal vector, the first order result is 

Sa^' = -^ 1 + hj\v 1 + E' 1 ). (C.9) 

(u) 

While this is still not equal to any version of ip directly, da) ' is, in fact, equal to Cl T, when 
the latter is evaluated using a uniform density slicing. Going to second order, we find 


s4 u) ' = -^2 + ^V 2 (u 2 + E' 2 ) + i ( -4C^Cf + 20 1 V 2 ui + + v() 


+ C h 


i v vi + v \ + + up* + B\i — S\i) 

1 


+ 


— 2(u|/i + v[ ) ( — ipi + —V / (Ei + vi)dr 


(C.10) 


Again, this variable is not equal to any version of '(/>, but it becomes exactly V’/, when evaluated 
using a comoving threading (v = Vy = 0). This is equivalent to saying that, by applying the 
same procedure to this quantity, one would obtain a gauge invariant quantity that is equal 
to C/ 2 '• This is not surprising, given the results of Refs. [11, 31, 35, 36], which found similar 
evolution equations for gauge invariants defined from the expansion scalar, 0 = 

We conclude our exposition of this appendix by noting that, even though the connec¬ 
tion between the intrinsic curvature and -0 is lost at second order, it is still possible to find a 
definition of ^ which closely matches the extrinsic curvature scalar for comoving hypersur¬ 
faces (i.e. space-like hypersurfaces which are normal to the velocity 4-vector), at least when 
evaluated in a particular gauge. The reason why the version of that resembles K is the one 
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arising from the determinant of g lJ can be explained by a relation between the determinant 
of the metric and the covariant divergence of a 4-vector. This is given by 

= d tl y! J + T ^u v = d tJ u IJ ' + u v d v log (V~g) , (C.ll) 

in which g = det [g^v]- Furthermore, it can be shown that det[( 7 * J ] is related to g by 

9 = 9oo (det^]) 1 , (C.12) 


and thus the previous relation becomes 

Vp/t/'' = + u v d v log (V~9oo) 


u u d u log 



(C.13) 


Choosing a comoving threading is equivalent to setting u l = 0 and in that case it is straight¬ 
forward to show that u° = y/—goo- This implies that the first two terms on the r.h.s. of 
Eq. (C.13) cancel, and one finds 

( V //«''')com = -u^log (j/detigvjj = 3^ (log a+ ipi) , (C.14) 


in which we substituted det \g l] ] by the definition of ipj. Equivalently, one has 


i [ K {u) ds 


= log a + ijji. 


(C.15) 


This shows V’/ to be the perturbation to the integrated extrinsic curvature of comoving 
hypersurfaces when written using a comoving threading. This result is valid at all orders 
and provides a clear interpretation to this perturbation derived from the determinant of the 
spatial part of the inverse metric. 
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